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Abstract 

We consider the three flavor Nambu-Jona-Lasinio model with the 't Hooft 
interaction incorporating the U(1)a anomaly. In order to set the coupling 
strength of the 't Hooft term, we employ the topological susceptibility x in- 
stead of the rj meson mass. The value for x ls taken from lattice simulations. 
We also calculate x a t finite temperature within the model. Comparing it 
with the lattice data, we extract information about the behavior of the U(1)a 
anomaly at finite temperature. We conclude that within the present frame- 
work, the effective restoration of the U(1)a symmetry does not necessarily 

take place even at high temperature where the chiral symmetry is restored. 
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I. INTRODUCTION 



The U(1)a anomaly of QCD plays an essential role in hadron physics. One of its most 
striking manifestations would probably be the rf meson mass. Since the U(1)a symmetry is 
broken not spontaneously but explicitly by the anomaly, rf cannot be regarded as a nearly 
massless Nambu-Goldstone boson like the other psuedoscalar mesons. In fact, the mass of 
rj is as large as the nucleon mass, i.e. m v i = 958 MeV. This is the so-called U(1)a problem. 

The topological susceptibility, \i is an essential quantity in considering the U(1)a prob- 
lem because it is related to through the Witten-Veneziano mass formula 



2iVf 2 2 2 

-pX = m n + m r,> ~ 2m Ki I 1 ) 

J TT 

where Nf — 3 is the number of the flavors and f n is the pion decay constant. This formula 
has been confirmed by calculating \ directly on the lattice J7|: The calculations give \ ~ 
(180 MeV) 4 , which is consistent with the value obtained by plunging into the formula (|]) 
experimental values of the pion decay constant and the meson masses. Thus, the topological 
susceptibility could tell us as much information about the U(1)a anomaly as does m v >. 

The tool we will employ in this work for the investigation of the U(1)a problem is the 
three flavor Nambu-Jona-Lasinio (NJL) model @-§| which can be used as an effective theory 
of QCD. The NJL Lagrangian we adopt here is 

C = C + U + A;, (2) 
£o = q(il^ - m)q } (3) 
C A = Gj2[(q\ a q) 2 + (qz l5 X a q) 2 }, (4) 

a=0 

£ 6 = —K [det q(l + lb )q + det q(l - j 5 )q] , (5) 
where the quark field q is a column vector in the color, flavor and spinor spaces, and A a is 



the Gell-Mann matrix in the flavor space with A = y2/3 diag(l, 1, 1). The determinants in 
Eq. (KI) are with respect to the flavor indices. 
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The free quark Lagrangian Cq contains the current quark mass term with m = 
diag(m u , mj, m B ), breaking the U(3)l <8> U(3)r symmetry explicitly Throughout this ar- 
ticle, we assume the exact isospin symmetry, i.e. m u = rrid- The term £4 generates the 
four-point couplings and is invariant under the U(3)l <8> U(3)r transformation. The six- 
point, determinant term £ 6 is what is called the 't Hooft interaction and breaks the U(1)a 
symmetry. This interaction simulates the U(1)a anomaly in our scheme, and the effective 
coupling constant K measures its strength. 

Let us review the status of parameter setting in the NJL model. The parameters to be 
fixed are, the current quark masses (m n = m^, m s ), the three- momentum cut-off (A), and 
the effective coupling constants (G and K) . The physical quantities usually taken as inputs 
are, m^, f n , m^, and m v >. The question we bring out here is as to m v r which has been used 
for the determination of K. As is well known, the NJL model lacks in confinement, and in 
fact, in this model 7/ decays into asymptotic qq states due to its large mass. Thus m v r in the 
NJL model is not a well-defined quantity. The alternative quantity for the determination 
of K we will use here is the topological susceptibility, which contains the information about 
the U(1)a anomaly and whose value has been given by lattice Monte-Carlo simulations, as 
mentioned above. One of the main purposes in the current work is to derive the expression 
for the topological susceptibility within the framework of the NJL model, and to fix the 
parameter K by means of % as an input. 

Recently, the behavior of the effect of the U(1)a anomaly at finite temperature has been 
discussed intensively [p|,pl]-[T^| . In particular, special attentions have been paid to whether 
or not the effective restoration of the U(1)a symmetry and the chiral phase transition occur 
simultaneously. This question is still controversial and is not settled yet. Here, we should 
clarify what we mean by 'the effective restoration of the U(1)a symmetry' It means that 
all U(1)a violating quantities vanish, i.e. that all order parameters of the U(1)a symmetry 
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vanish.[] The possibility that rf would degenerate with the other pseudoscalar mesons has a 



great deal of significance upon the experimental view in relativistic heavy ion collisions [|T0 
In the NJL model, the definition we gave above is equivalent to K getting to zero, since 
finite K makes the U(1)a symmetry breaking manifest in the NJL lagrangian (JJ). Since the 
origin from which a finite value of K arises is the presence of instantons in the physical state, 
the effective restoration of the U(1)a symmetry is expected owing to the naive argument 
that the instanton density will be suppressed at sufficiently high temperature. Thus when 
the magnitude of K becomes smaller, we will call it 'the effective restoration of the U(1)a 
symmetry' 

The temperature dependence of K in the NJL model, which indicates nothing but the 
temperature dependence of the U(1)a anomaly, has been put by hand and not gone beyond 
phenomenology [|J. This is both because experimental data for at finite temperature, 
which are necessary for determination of K, are not available, and because rf becomes 
unbound completely in the model soon after we raise the temperature from zero. Since the 
topological susceptibility has been calculated at finite temperature on the lattice, we will 
be able to determine the temperature dependence of K using that data, and obtain some 
knowledge about the effective restoration of the U(1)a symmetry. 

The paper is organized as follows. In Sec. J|, we will derive the expression for x m the 



NJL model. Sec. Ill is devoted to parameter setting of the model and numerical calculations 



of the physical quantities. A summary is given in Sec. IV 



II. TOPOLOGICAL SUSCEPTIBILITY IN THE NJL MODEL 

In this section, we calculate the topological susceptibility, Xi within the framework of 
the NJL model. 



x The authors thank the unknown referee for his suggestion on the definite meaning of the U(1)a 
symmetry restoration. 
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The first task is to know a general expression of x- We recapitulate here the definition 
of x- We begin with the QCD Lagrangian density 

£qcd = ~\f; v F^ + q(ij^ - m)q + 6Q, (6) 

where F* is the gluon field strength tensor, = + igA^ is the covariant derivative 
with being the gluon field, 9 is the QCD vacuum angle, and Q is the topological charge 
density defined by 

With this Lagrangian density, the vacuum energy density e is written as 

e~ £VT = J VA^DqDqeS dixC ^ = Z (8) 

in a path integral form, where V and T are the space and time volumes, respectively. The 
topological susceptibility x is defined as a second derivative of e with respect to 9 at 9 = 0, 

= ^1 
X ~ d9 2 

where T stands for the time ordering operator, and the subscript 'connected' means to pick 
out the diagrammatically connected contributions. Thus, in order to calculate x i n the NJL 
model, it is necessary to find a correspondent to Q(x) in the model. For that purpose, we 
consider the four- divergence of the U(1)a current, Js M = <?7 M 75<?. In QCD, one has 

&>J 5li = 2NfQ(x) + 2tqm l5 q, (10) 

which does not vanish due to the anomaly. 

On the other hand, in the NJL model (Eqs. (@)~ (§[)), we find ||] 

d» J 5fl = AN { Klm det $ + 2iqm-f 5 q, (11) 

where 

= -75)^, (12) 



/ d 4 x(0|TQ(x)Q(0)|0) 

connected j (9) 



and i,j denote the flavor indices. Comparing these two expressions, we find that 

Q(x) = 2inmdet <P 

= -iK [det * - (det $)*] (13) 

in the NJL model. 

With the definition of r± = 1 ± 75, we can write 

det <p = y^^q i r_q a )[^r.q b ){q h r.q c ), 
(det #)• = ^6^6 /mn (g / r + g d )(g m r + g e )(g-„r + g / ) (14) 

so that 

connected 

K 2 r 

— / j4 abc ijk def Imn 

~ (3!) 2 7 

x (0| T {(gjr_g a ) (q^T^qb) (q k r_q c ) (x) (g/T_g d ) (q m r^q e ) (q n r^q f ) (0) 
-(qir-q a )(qjr_q b )(q k r_q c )(x)(qir + q d )(q m r + q e )(q n r + q f )(0) 
-(qir + q a )(qjr + q b )(q k r + q c )(x)(qir _q d )(q m r -q e )(q n r -<?/)(0) 
+ (g i r + g a )(g j r + g 6 )(g fc r + g c )(x)(gir + g d )(g m r + g e )(g:„r + g / )(0) } |0) connC cted- (15) 

Now we must evaluate these four matrix elements. For the time being, we pick up one term 
out of the four. Following Wick's theorem, we take full contraction in terms of the propagator 
S(x, x') that has been constructed in the self-consistent gap equation 0-§J. Although several 
ways of contraction are possible, there exists only one that is the leading order in expansion 
in terms of inverse powers of the number of colors, i.e. 1/N C . The situation is demonstrated 
in Fig. ^ by means of the finite range representation. 

The diagrams (a) and (b) in Fig. [I] contain iV c 4 coming from traces over color. The 
diagram (c) contains iV c 5 and is the leading order in 1/N C expansion. Notice that (b) is the 
exchange term for (c) and is lowered by 1/N C as compared with (c). Since the gap equation 
for the constituent quark masses and the dispersion equations for the meson masses are 
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derived up to the leading order of the large N c expansion ||, we should take only the 
contribution of (c) for the consistent treatment. 

Taking account of the four terms in Eq. flT5|), we obtain the following expression for the 
lowest order of the diagrammatical expansion, 

X (1 ° wcst) = -^(-9)e abc e^ k e def e lmn 4y d A xN c ti [S dl (x) l5 S al (x) l5 ^ 

xN c Hr [S bj (0)\ tr [5^(0)] tr [S em (0)] tr [S fn (0)} , (16) 

where the full propagator in the Euclidean space is given as 

S-.(x)-8- f d4p t + m * c -ivx (17) 

and m* denotes the constituent quark mass. 

The object in the curly brackets in Eq. (|lfj) corresponds to the one-loop part connecting 
the points x and in Fig. |TJ(c). This is basically the one-loop proper polarization insertion 
nfj {k 2 ) @~|6[] with k 2 = although the trace over flavor is not taken in this case. Especially 
in the case of a — i and d = I, which is the condition for giving non-zero contribution to 
^.(lowest) ^ ue ^ ^ ne £_t ensor i n Eq. ([TB]), the object can be identified with an element of 3-, 
8-, 0-channel polarizations with k 2 = 0; for instance, 0-0 channel polarization is 

d 4 p 



n P oo(k 2 = 0) = tr flavol | diag (-N c J 



(27T) 



tr[z 75 5 u (p) ?7 5^ u (p)] 



d P . r- od/ \ • od/ <\ A r f d P 



N r / — 4rtr 



t~? 5 S d (p)t l5 S d (p) , -N c / ^tr [i l5 S s (p)i l5 S s (p)} . (18) 



(2tt) 4 L '° ^' 10 ^'i' L J (2vr) 4 ' 
Actually, there exist other diagrams that are of the same order in 1/N C expansion as the 
diagram in Fig. [l|(c). They are shown in Fig. 0. They are of the same order as Fig. |I|(c) 
because, while each four-point vertex is of 0(N C ^ 1 ) 0, it is compensated by a factor iV c 
coming from its neighboring loop. We have to include all these diagrams for consistency of 
the 1/N C expansion. We will call their sum x^ rmg \ f° r they are regarded as the ring diagrams 
to be resummed in the mean field approximation. 

Of course, the sum of these ring diagrams with the one-loop diagram included can be 
interpreted as a propagation of a certain meson. Note that the momentum of the propagating 



particle is zero; k 2 = 0. This is just the reflection of the fact that x is the quantity of the 
zero frequency mode of the Fourier transform of (TQ(x)Q(O)), namely, 

x = f d i xe-^(TQ(x)Q(0)) . (19) 

Now we calculate ^( rin s) following the Feynman rules of the NJL Lagrangian. We note 
that the diagrams in Fig. § are obtained by replacing the one-loop polarization part in Fig. 
gKc) with the ring diagrams. Correspondingly we can obtain ^( rmg ) by replacing the curly 
bracketed part in Eq. (|16|) with the sum of the ring diagrams. We first perform the sum- 
mation over flavor indices. After that, it can be shown that the expressions corresponding 
to each edge point of a diagram in Fig. ^| are brought together into a matrix form and are 
arranged to the linear combination of Xs and Ao matrices. In other words, the 8- and 0- 
channel vertices have been assigned to each edge point of each ring diagram in Fig. |2|. Then 
following the Feynman rules, we can construct the ring diagrams by linking 8-, 0-channel 
polarizations UH(k 2 = 0), ITfo(£; 2 = 0) = n^(fc 2 = 0), hQ£; 2 = 0), and 8-, 0-channel vertices 
Kg% \ = Kq$ \ Kqq , in all possible ways. The result is 

x2K(l- 2TLK) ~ 



where 



x j_L trS ~ (trS . _ trS «) g j + J_ trS » (2tr5 - + trS u } gj 1 (20) 



n.f") (2!) 
V k$ I V n 08 n 00 / 



are 2x2 matrices, and 



= U^(k 2 = 0), tiS 1 = tiS\x = 0). (22) 

For example, the first term in the first curly brackets of Eq. (|20"D produces contributions 
from the diagrams whose one edge point is an 8-channel vertex. And, again for example, 
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the second term in the second curly brackets of Eq. (20) produces contributions from the 
diagrams the other edge point of which is a 0-channel vertex. The 8-8 and 0-0 channel 
diagrams can be interpreted as the propagations of r}$ and rjo mesons respectively. We see 
that in addition to the 7/g and r/o propagations, there occur 8-0 and 0-8 mixing channel 
diagrams. 

Another comment is in order. If trS u = trS s , that is, the SU(3)v symmetry is exact, the 
8-channel vertices in Eq. (|20| ) vanish. In this case, 8- and mixing channel polarizations ITgg, 
llgo = flgg as well as 8- and mixing channel vertices i<Qg , = all vanish, so that 
only the ring diagrams constructed by those of the 0-channel, UqI and , contribute to 
^( rm s) ? which are interpreted as the propagation of the pure rj state. 

Finally, combining Eqs. ([16]) and (|20|), we arrive at the expression for the topological 
susceptibility, 



Let us give one more comment. In general, a two-point correlation function of gauge 
invariant operators may be decomposed into the sum over multi-particle intermediate states 
by inserting a complete set between two operators. It is known that the dominant contri- 
butions of the leading order in 1/N C expansion are those of one-particle intermediate states. 



the momentum of the intermediate particle is zero, which is the case of x> the contribution 
of the rjo propagating state is the only leading term in 1/N C expansion. These statements are 
quite consistent with our specific model calculation respecting large N c expansion (although 
r/s and mixing channels besides rj propagate in the intermediate states in our model due to 
the explicit SU(3)v symmetry breaking). 




(23) 



Moreover Witten fJTJ| and Veneziano have derived their formula by assuming that, when 
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III. NUMERICAL CALCULATION 



Now that we have obtained the expression for the topological susceptibility, we proceed 
to numerical calculations. In Sec. Hi A , we set the parameters at zero temperature by 
employing x- With the determined parameters, we calculate physical quantities. In Sec. 
Ill Bl we discuss the temperature dependence of x an d the six-point coupling constant K. 



A. Parameter setting at zero temperature 



The parameters to be set in the NJL model are; 

current quark masses 
three-momentum cut-off 
four-point coupling constant 
six-point coupling constant 



A 
G 
K. 



As for m u = m&, we set them to be m u = m& = 5.5 MeV following Ref. @]. 
To set the other four parameters, we use the following quantities as inputs, 

= 138 MeV 

U = 93 MeV 

m K = 495.7 MeV 

X = (175±5MeV) 4 . 

The fourth quantity we use here in place of m v r = 957.5 MeV is, as mentioned in the 
introduction, the topological susceptibility x- The numerical value of x is taken from Ref. 
0, in which x is calculated in the quenched approximation. 

Initially, however, we will calculate with the parameters determined by using m v t as 
input in order to check consistency of and x i n the NJL model. Parameter setting with 
rrirf has been performed in Ref. M, and the results are 



m s = 135.7 MeV, A = 631.4 MeV, GA = 1.835, KA = 9.29. 
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The physical quantities calculated from these parameters are summarized in the first column 
of Table |. We first check the Witten-Veneziano mass formula (|IJ) within the NJL model. 
The computed values of \ and m v with those parameters of Ref. are 
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Xnjl = (166 MeV) 
m v = 487 MeV, 

so that the ratio of LHS to RHS in Eq. (fj) turns out to be 

n( 2 ™f 9 2 , = 0.81. (24) 
/2( m 2 +m 2, -2m 2 K ) 

On the other hand, the ratio of Xnjl to Xhat, which means how much the conventional 
parameters determined with m v r reproduces the lattice data of Xi is 

^ = 0.80. (25) 

XLat 

^From the above two ratios, we can say that as a whole, and % are reproduced well 
simultaneously in the NJL model. 

Now we consider parameter setting with x used. The topological susceptibility, x might 
be a more suitable quantity for parameter setting than m v i for the following two reasons: 

• Since rj' decays into asymptotic qq state due to lack of confinement in the NJL model, 
m n i may be a less reliable quantity, while x is f ree from such a shortcoming of the NJL 
model. 

• The value of Xi (175 MeV) 4 is small enough compared with the cut-off A ~ 600 MeV. 
Thus the NJL model is expected to describe x we U- 

The parameters obtained by using x — (175 MeV) 4 are 

m s = 135.7 MeV, A = 631.4 MeV, GA 2 = 1.765, KA 5 = 11.32. 

We note that KA 5 becomes larger than the case of using m v i as an input, which implies 
that the binding of rj' is loosened. (The 't Hooft interaction loosens the binding of mesons, 
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that is, induces a repulsive force between quarks. This can be seen from the very fact 
that rf becomes massive due to the interaction.) Physical quantities calculated with these 
parameters are shown in the second column of Table |. The solution for m v i in the mean 
field approach does not exist, that is, rj is not bound any more. We see that m v is improved 
slightly. Although rj no longer exists in the NJL model, we could infer its mass by utilizing 
the Witten-Veneziano mass formula ([[]); m v ' = 942 MeV is obtained. 

We close this subsection by referring to the study due to Takizawa, Nemoto, and Oka 
||, in which the parameters, especially the six-point coupling constant K, are determined 
in a different approach. They examined the radiative decays of an rj meson such as rj — > 27, 
f] — > 7i~r l ~ and rj — > 7r°77, and obtained rather strong six-point coupling constant KA 5 , 
namely, four times as large as that determined by using m v i. Although we cannot compare 
our parameters directly with theirs due to different cut-off schemes (their scheme is the 
four-momentum cut-off), it is not probable that our result is compatible with theirs. Still, 
we believe that our approach is rather straightforward to probe the U(1)a anomaly. 

B. Behavior of K at finite temperature 

In this subsection, we discuss the temperature dependence of K, comparing the NJL 
calculation of x with the lattice data. 

The lattice data for the topological susceptibility [[F] are shown in Fig. [3] with error 
bars. The T c in the figure denotes the temperature of the chiral phase transition. Although 
T c = 260 MeV in the original Ref. 0, we have rescaled it to 150 MeV. We should notice 
that the lattice data are computed only up to T = 1.4T C . Unfortunately, the lattice data are 
absent at high temperatures. At any rate, the data show that x drops rapidly around T c . 

One comment should be noted. The fact that x drops near T c does not always mean 
the effective restoration of the U(1)a symmetry at T c . This can be seen by returning to the 
Witten-Veneziano mass formula ([!]), 

2N iX = f*(m 2 v + m 2 v ,-2m 2 K ). (V) 
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We realize that the pion decay constant /„-, which is associated with the spontaneous chiral 
symmetry breaking (SU(3)l <8> SU(3)r — > SU(3)v), has entered the formula. Since /„. would 
become zero along with the restoration of the chiral symmetry, \ is also expected to become 
zero around T c . In this sense, the lattice data which show the dropping of x & t T c is 
what should be expected from the formula, and rather, we could consider that the data 
confirm the validity of the Witten-Veneziano mass formula. Thus the dropping of x i n the 
lattice data should be attributed to the restoration of the chiral symmetry, and does not 
always indicate the effective restoration of the U(1)a symmetry. It is worth noting that 
this behavior results from large N c expansion. In Ref. [[U| it was pointed out that n-point 
correlation functions (n < Nf = 3) cannot detect any effect of the U(1)a anomaly in the 
chiral symmetric phase. One might have thought that the dropping of x near T c would be 
regarded as x's insensitivity to the U(1)a anomaly. However, it is not the case because x is 
not a U(1)a singlet quantity. In fact, it contains contributions carrying the U(1)a charge 2, 
and -2. Thus x is an appropriate quantity to observe the fate of the U(1)a symmetry even 
in the chiral symmetric phase, in principle, even though careful attention should be paid in 
order to infer correct meanings. 

In this respect, the discussions of Schaffher-Bielich (JTSj is obscure; he discusses under 



the assumption that x an d the U(1)a anomaly are equivalent to each other and that the 
dropping of % & t T c immediately means the effective restoration of the U(1)a symmetry. His 
assumption is considered as correct only if the dropping rate of x is much faster than that 
of /„-. To judge the validity for this prevailing assumption is what we pursue in the present 
work. In fact, as discussed below, our result reveals that the assumption has no convincing 
reliability, at least, within the framework of the NJL model. 

Now we consider the temperature dependence of x i n the NJL model. Among the four 
parameters (m s , A, G, K), we might reasonably fix m s and A at the values determined at 
zero temperature. In general, however, we should take account of temperature dependences 
of the coupling constants GA 2 and KA 5 . As for GA 2 , it would be hard or almost hopeless to 
get information about the temperature dependence even in some phenomenological sense. 
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Here, we make an assumption that GA 2 does not depend on temperature. This might be 
partially justified by the fact that even if GA 2 is constant, the NJL model restores the chiral 
symmetry as a consequence of its own dynamics. 

We now pay attention to the behavior of KA 5 , which indicates the temperature depen- 
dence of the U(1)a anomaly. For the first case, we treat KA 5 as a constant parameter and fix 
it at the values at zero temperature. We will call this prescription CASE A. The calculated 
temperature dependence of x is shown in Fig. || (the solid line). We see that x i n the NJL 
model drops near T c as the dynamical consequence and reproduces the lattice data up to 
1.4T C considerably well. This result means that the U(1)a symmetry is not restored at least 
up to 1.4T C , and we conclude that the effective restoration of the U(1)a symmetry does not 
coincide with the chiral phase transition. 

At high temperatures, of course, we cannot judge whether or not the U(1)a symmetry is 
restored, since we lack the lattice data in those temperatures. If we believe that the instanton 
density is suppressed exponentially at high temperatures as is expected by the Pisarski- 
Yaffe factor 0, and the correlation of the topological charges, i.e., x is & lso suppressed 
exponentially, the fitted line for the lattice data with a Fermi function in Fig. |3| (the dashed 
line) could be considered as reasonable behavior. We notice here the deviations of the NJL 
calculation (the solid line) from the fitted line at high temperatures. As the CASE B, we 
let KA 5 have the temperature dependence such that it reproduces the fitted line of x- The 
calculated temperature dependence of KA 5 for this case is shown in Fig. [|. 

We notice that there are two lumps around T c . It would be senseless to take them 
seriously since we have ignored the temperature dependence of GA 2 that should have been 
taken into account in principle. Rather, we should note that the U(1)a symmetry is restored 
at high temperatures as is expected from the starting assumption that the instanton density 
is suppressed at those temperatures; the consistency is maintained in the NJL model. 

We now calculate the constituent quark masses, the meson masses, and the pion decay 
constant in our CASE A and CASE B. The results are shown in Figs. The qualitative 

features are almost the same as those of the CASE I by Hatsuda and Kunihiro @]. 
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Finally, we give the temperature dependence of m v i in Figs. ^ and [14] that could be 
obtained by utilizing the Witten-Veneziano mass formula flip. The rj mass goes to infinity 
at around T c = 200 MeV in either CASE A or B. This is because f w gets to zero at that 
temperature. We have removed those infinities above the temperature at which f v vanishes 
because it is considered that our approximation scheme is broken down there. 



IV. SUMMARY 

We have derived the expression for the topological susceptibility, Xi i n the NJL model 
within the same approximation as for the constituent quark masses and the meson masses, 
namely in the leading order of large N c expansion. At zero temperature, we have performed 
parameter setting by employing x i n place of m v > , and have seen that the obtained parameters 
do not allow the bound state of rj'. At finite temperature, we have calculated the temperature 
dependence of Xi an d have found that the lattice data up to 1.4T C are reproduced with 
a constant six-point coupling constant K. This means that the U(1)a symmetry is not 
restored up to 1.4T C , and we are led to the conclusion that within the present framework, the 
effective restoration of the U(1)a symmetry and the chiral phase transition do not necessarily 
occur simultaneously even though the rapid dropping of x around the chiral transition, 
observed in the lattice simulation, seemingly suggest the simultaneous restoration. At high 
temperatures, we cannot state anything definitely because of the absence of the lattice data. 
We have shown, however, that if x ls suppressed exponentially, the U(1)a symmetry is 
allowed to be restored at high temperatures. 

The topological susceptibility is an interesting quantity because it is related to the mass 
of rf through the Witten-Veneziano mass formula. At zero temperature, we have seen that 
the formula is satisfied numerically in the NJL model. At finite temperature, by utilizing 
the formula, we have obtained knowledge as to the temperature dependence of m v i. In NJL 
model, r)' is far from a stable particle even if it exists. Therefore we would say that the 
parameter setting by using my is somewhat obscure. Our approach proposed in this article 
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is, on the other hand, not affected by questionable quantities such as m n i and as a result it 
takes the advantage to extract reliable results even when the bound state of rf cannot be 
available in the NJL model. 
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FIGURES 

FIG. 1. Three of various contracting ways are shown with their corresponding diagrams. 

FIG. 2. Other leading order diagrams. 

FIG. 3. The lattice data are plotted with error bars. We have fitted them with a Fermi function 
(the dashed line). The solid line denotes the result of the NJL model with constant KA 5 . 

FIG. 4. Temperature dependence of KA 5 . 

FIG. 5. The constituent quark masses for CASE A. 

FIG. 6. The constituent quark masses for CASE B. 

FIG. 7. The meson masses for CASE A. 

FIG. 8. The meson masses for CASE B. 

FIG. 9. The Kaon mass for CASE A. 

FIG. 10. The Kaon mass for CASE B. 

FIG. 11. The pion decay constant for CASE A. 

FIG. 12. The pion decay constant for CASE B. 

FIG. 13. The rf mass for CASE A. 
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FIG. 14. The rf mass for CASE B. 
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TABLES 



Hatsuda and Kunihiro j|] Ours Experimental/empirical values 

m u * (MeV) 335 337 336 

m s * (MeV) 527 523 540 

m ri (MeV) 487 505 549 

m v , (MeV) (958) None [942] 958 

X 1/4 (MeV) 166 (175) 175 

V -21° -16.7° -20° 

TABLE I. Calculated physical quantities. Comparison of our results with Hatsuda-Kunihiro 
and experimental data. In parentheses are the values used as inputs. The value of m„/ shown in 
the square brackets is inferred by the Witten-Veneziano mass formula. 
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